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Abstract 

The turbulent transport of main ion and trace impurities in a tokamak device in the presence 
of steep electron density gradients has been studied. The parameters are chosen for trapped 
electron (TE) mode turbulence, driven primarily by steep electron density gradients relevant 
to H-mode physics, but with a transition to temperature gradient driven turbulence as the 
density gradient flattens jj] . Results obtained through non-linear (NL) and quasilinear (QL) 
gyrokinetic simulations using the GENE code [2)01 compared with results obtained from 
a fluid model. Main ion and impurity transport is studied by examining the balance of 
convective and diffusive transport, as quantified by the density gradient corresponding to 
zero particle flux (peaking factor). Scalings are obtained for the impurity peaking with the 
background electron density gradient and the impurity charge number. It is shown that the 
impurity peaking factor is weakly dependent on impurity charge and significantly smaller 
than the driving electron density gradient. 

The compatibility between a reactor-grade plasma and the material walls surrounding the plasma 
is one of the main challenges facing a magnetic fusion device. The presence of very low levels 
of high Z impurities in the core plasma may lead to unacceptable levels of radiation losses and 
fuel dilution. Also low Z impurities, in the form of Beryllium or Helium-ash, may result in 
fuel dilution that severely limits the attainable fusion power Consequently, the transport 
properties of impurities is a high priority issue in present experimental and theoretical fusion 
plasma research. This is emphasised by the the new ITER-like wall experiment in JET ^] , where 
a Beryllium-clad first wall in the main chamber, combined with carbon and tungsten tiles in the 
divertor, will be tested for the first time. 

The transport of main fuel as well as impurities in the core region of tokamaks is expected 
to be dominated by turbulence driven by Ion Temperature Gradient (ITG) modes and Trapped 
Electron (TE) modes. The main drives for the ITG/TE mode instabilities are gradients of 
temperature and density combined with unfavourable magnetic curvature. Most of the theoret- 
ical studies of turbulent particle transport have been devoted to temperature gradient driven 
ITG and TE modes, using both fluid, quasilinear (QL) and nonlinear (NL) gyrokinetic mod- 
els Much less effort has been devoted to particle transport in regions with steep density 
gradients. The density gradient provides a drive for TE modes which may dominate the tem- 
perature gradient drive for plasma profiles with i?/L„_. > R/Lt^. This may occur in connection 
with the formation of transport barriers, like the high confinement mode edge pedestal, in fusion 
plasmas. 

In the present letter, the turbulent transport of main ion and trace impurities in tokamaks 
is investigated through nonlinear (NL) gyrokinetic simulations using the GENE code. The main 
part considers collisionless TE modes driven by density gradients but with a transition to tem- 
perature gradient driven TE modes as the density profiles flattens. The impurity density gradient 
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for zero impurity flux is calculated for varying background electron density gradient drive and 
for a range of impurity species. This study complements recent studies 2^ 27 1 on temperature 



gradient driven TE and ITG mode impurity transport. The results are compared with QL ki- 
netic simulations and a computationally efflcient multi fluid model, suitable for use in predictive 
transport simulations. Of particular interest is the sign of the impurity convective flux and the 
degree of impurity peaking in the presence of strong background electron density gradients. 

The models used have been described in detail elsewhere, see [2^ and references therein, only 
a brief summary is given here. The NL and QL GENE simulations were performed in a flux tube 
geometry, in a low /3 (/3 = 10"'') s-a equihbrium [I. I28l - l30l | . 

In order to ensure that the resolution was adequate, the resolution was varied separately for 
the perpendicular, parallel and velocity space coordinates, and the effects of this on the mode 
structure, k± spectra and flux levels were investigated. The resolution was then set sufficiently 
high for the effects on these indicators to have converged. For a typical NL simulation for main 
ions, fully kinetic electrons, and one trace species, a resolution of n^: x rij, x = 96 x 96 x 24 
grid points in real space and of n„ x = 48 x 12 in velocity space was chosen. For QL GENE 
simulations the box size was set to n^; x x = 8 x 1 x 24 and riy x = 64 x 12 respectively. 
The impurities were included self-consistently as a third species in the simulations, with the trace 
impurity particle density n^/ng = 10^^ in order to ensure that they have a negligible effect on 
the turbulence. 

For the fluid simulations, the Weiland multi-fluid model [3l[ is used to derive the main ion, 
impurity, and trapped electron density response from the corresponding fluid equations in the 
collisionless and electrostatic limit. The free electrons are assumed to be Boltzmann distributed. 
The equations are closed by the assumption of quasineutrality. An eigenvalue equation for TE and 
ITG modes is obtained in the presence of impurities. The eigenvalue equation is solved for general 
mode width 31|. Alternatively, a strongly ballooning eigenfunction with k?, = (Sq^R^) can 



be used for magnetic shear s ~ 1 [3^. The eigenvalue equation is then reduced to a system of 
algebraic equations that is solved numerically. 

The main ion and impurity particle fluxes can then be written as: 

r, = {SujVExB) = -rijPsCs ■ (1) 

The angled brackets imply a time and space average over all unstable modes. Performing this 
averaging for a fixed length scale keps of the turbulence, the particle fiux can be written: 



-D,^ + Dt,^+RV,,. (2) 



The first term in equation Q corresponds to diffusion, the second to the thcrmodiffusion and 
the third to the convective velocity (pinch), where 1/ Ln. ~ —Vrij /rij, nj is the density of species 
j and R is the major radius of the tokamak The pinch contains contributions from curvature and 
parallel compression effects. These have been described in detail in previous work 20-22, 2(|. 
For trace impurities, equation Q can be uniquely written RFz/nz = DzR/Ln^ + RVz, where 
Dz is the impurity diffusion coefficient and Vz is the impurity convective velocity. 

The zero- flux impurity density gradient (peaking factor) is defined as PFz = —RVz/Dz for 
the value of the impurity density gradient that gives zero impurity fiux. Solving the linearised 
equation ^ for R/Lm with Fz = yields the interpretation of PFz as the gradient of zero 
impurity flux, and it quantifies the balance between convective and diffusive impurity transport. 

The main parameters used in the simulations are summarised in table [T] The parameters 
where chosen to represent an arbitrary tokamak geometry at about mid radius, and do not 
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represent any one particular experiment. A a moderately steep electron temperature gradient 
{R/Lt, = 5.0) together with a flatter ion temperature gradient {R/Lxi^z = 2.0) were used to 
promote TE mode dominated dynamics. Following the background density gradient for the 
base scenario was set higher than the temperature gradient, to ensure density gradient driven 
dynamics. In order to preserve quasineutrality Vrig — Vn^ was used. 

First, the main ion particle flux (Fp) is studied. Time averaged fluxes are calculated from 
time series of NL GENE data after convergence, as illustrated in figure [Ta] The scalings of Fp 
with the electron density gradient obtained from NL GENE and fluid simulations are shown in 
figure llbl The large transport found in NL GENE simulations is an indication of the stiffness 
of the gyrokinetic model, and is often seen in fixed-gradient simulations of turbulence. The fluid 
model shows a similar scaling of the main ion flux, but the transport is smaller for i?/i„_, > 3.0. 
The main ion density gradient corresponding to zero ion flux {PFp) can be found by similar 
means to that of PFz, however, since the trace approximation is not valid for the main ions, 
the zero-flux gradient has to be found explicitly by varying Vnp until the condition Fp ~ is 
satisfied. The NL GENE results presented in figure [Tbl indicate that for the present parameters, 
lower density gradients only result in Fp — > 0. The the fluid model gives a small outward flux in 
the limit R/Ln^ — > 0. Neither model results in flux reversal for TE mode driven turbulence. 

Next, the scaling of the impurity transport with the background density gradient (i?/L„^) 
is investigated. The results for the impurity peaking factor are shown in figure I2al We note 
that the impurity peaking saturates with PFz ~ 2.0 for large values of the electron density 
gradient. The QL results tend to consistently overestimate the peaking factors, while the fiuid 
model gives results that are somewhat below the NL GENE results for the steeper gradients. 
The fluid results show a considerably less dramatic dependency of the peaking factor than the 
gyrokinetic results, both of which show a strong decrease in PFz as the electron density profiles 
flatten. This is observed for all values of the impurity charge number. As the background density 
profile becomes more peaked, a corresponding increase in impurity transport is expected. This 
is illustrated in figure [2bl where scalings, obtained from NL GENE simulations, of the diffusivity 
{Dz) and convective velocity (RVz) with R/Ln^ are shown. Although Dz and RVz strongly 
increase with R/Ln^, the impurity peaking {PFz = —RVz/Dz) is only weakly sensitive to the 
electron density gradient. For R/Ln^ ^ 2.0 the impurity peaking factor is not well defined, since 
both Dz and RVz go to zero. The corresponding linear eigenvalues are displayed in figurej^cl The 
fluid and gyrokinetic results are in qualitative agreement, showing an growthrate that increases 
uniformly with _R/L„_,. The results indicate a smooth transition from density gradient driven to 
temperature gradient driven TE mode turbulence, which dominates for R/Ln^ ^ R/Lt^ [ij. 

The scaling of the impurity peaking factor with impurity charge (Z), with _R/i„_, as a pa- 
rameter, is illustrated in figure [31 The models show only a very weak scaling, with PFz falling 
toward saturation for higher Z. The results are similar to those for the temperature gradient 
driven TE mode reported in Notably, the QL GENE simulations overestimate the peaking 
factors, whereas the fluid results are lower than the peaking factors obtained from NL GENE 
simulations. The trend observed for low Z impurities is reversed compared to trends reported in 
e.g. [26[ for ITG mode driven impurity transport. The qualitative difference can be understood 
from the Z-dcpendent thermodiffusion in equation which is outward for ITG modes and 
inward for TE modes. 

In summary, the turbulent transport of main ion and trace impurities in regions of steep 
density gradients has been investigated through nonlinear (NL) gyrokinetic simulations using the 
GENE code. The main part has considered coUisionless TE modes driven by density gradients 
but with a transition to temperature gradient driven TE modes as the density profiles flattens. 
The results for the impurity density gradient of zero particle flux (peaking factor) have been 
compared with QL kinetic simulations and a reduced and computationally efficient multi-fiuid 
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model, suitable for use in predictive transport simulations. 

For the parameters studied, qualitative agreement between gyrokinetic and fluid results has 
been obtained for the scaling of the impurity peaking factor with both the background density 
gradient and the impurity charge. In the region of steep electron density gradients, it was shown 
that the impurity peaking factor saturates at values significantly smaller than the driving electron 
density gradient. It was noted that for the chosen length scales, the QL GENE results generally 
overestimate the peaking factor, whereas the fluid results are close to or lower than the NL GENE 
results. The scaling of the peaking factor with impurity charge was observed to be weak, with a 
slight increase in the impurity peaking factor observed in the gyrokinetic results for low impurity 
charge numbers. 
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Table 1: Parameters used in the gyrokinetic simulations, ^ denotes scan parameters 
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(a) time series and time averages of the main ion fiux (Fp) from NL GENE simulations 
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(b) main ion flux (Fp) dependence on the background density gradient {R/Ln,)- 

Figure 1: Main ion flux (Fp) dependence on the background electron density gradient 
{—RVne/rie = R/Ln^)- NL GENE and fluid data with protons as main ions. Parameters 
are q = 1.4, s = 0.8, e = r/R = 0.14, R/Lt„z = 2.0, R/Lt, = 5.0, and r ^ Te/T; 1.0. The 
fluid data was obtained for kgps = 0.2. The fluxes are normalised to VT^inepf / R^ ■ The error 
bars indicate an estimated uncertainty of one standard deviation. 
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(a) dependence of the impurity peaking factor (PFz) on the background density 
gradient 
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(b) dependence of the impurity diffusivity and convective velocity {Dz and RVz) on 
the background density gradient 
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(c) scaling of real frequency (cj,.) and growthrate (7) with the background density 
gradient 

Figure 2: Scalings of the impurity peaking factor {PFz — ~RVz / Dz) with the background 
electron density gradient {R/LnJ, with parameters as in figure [TJ QL and fluid data have been 
acquired using kgps = 0.2. Figure [2b] shows the diffusivities and pinches corresponding to the 
NL GENE impurity peaking factors (PFz) in figure!^ Dz and RVz are normalised to VT,iPi/R- 
The eigenvalues in figure \7c\ are from fluid and GENE simulations, and are normalised to Cg/R. 
The error bars indicate an estimated uncertainty of one standard deviation. 
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Figure 3: Scaling of the impurity peaking factor {PFz = —RVz/Dz) with impurity charge Z, 
with parameters as in figure [TJ kgps — 0.2 was used in the QL and fluid simulations. The error 
bars indicate an estimated uncertainty of one standard deviation. 
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